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1. Introduction 

Two-particle correlations provide information about the space-time structure and dynamics of 
the emitting source [jl|, |2[ |3|, ^ Considering the correlations that occur in relativistic heavy-ion 
collisions one usually assumes that: (i) the particles are emitted independently (or the source is 
completely chaotic), and (ii) finite multiplicity corrections can be neglected. Then the correlations 
reflect a) the effects from symmetrization (anti-symmetrization) of the wave function and b) the 
effects that are generated by the final-state interactions of the detected particles between themselves 
and with the source. 

The nominal quantity expressing the conelation function in terms of experimental distributions 
[|]is 

where A (k) = E(PN/(Pk and P2{KM) = EaEb d^N/id^kad^h) are single- and two-particle 



cross-sections and £",• = (o{ki) = ^ rrfi + k? is free particle energy, i = a,b. 

From the very first experiments while analyzing experimental data the people wanted to sepa- 
rate the final state interactions (FSI) from that which is due just to the reaction zone dynamics, as 
it was thought. At first sight, the FSI can be regarded as a contamination of "pure" particle cor- 
relations. However, it should be noted that the FSI depend on the structure of the emitting source 
and thus provide information about source dynamics as well ||^. Appreciable understanding of the 
particle-particle and source-particle FSI was achieved last years g Ig |T|, [11], 0, Mean- 



while, there are several problems connected to the distortion of particle-particle FSI which are not 
clear enough so far. In the present paper we shall focus on two of them. First one is an old problem 
of the strong final state interactions. The second one occur in the last years due to the growth of 
the energy of colliding nuclei and consequently is due to impressive increasing of the secondary 
particle multiphcities (1000 secondary pions at SPS and 2000-3000 at RHIC). 

Strong FSI is a complex problem because on the fundamental level it concerns the non- 
perturbative QCD. Solution of the problem in this aspect is beyond the scope of the paper, we 
solve the problem in the model approach just to have reasonable estimations of the influence of 
strong two-particle interaction (pions hard core repulsion) on the correlation function. In different 
approaches it was discussed in Refs. [14, 15, 16, [I^, 18] (see also [^, Section 3.2). 



Actually, the release of pions from a small spatial region in any case is sensitive to hard core 
repulsion of registered particles. Technically this means the following: i. Particles created in 
finite volume can be represented by localized states, for the sake of simplicity let us say by cen- 
tered Gaussian packets; ii. To obtain correlations we integrate over all positions of two Gaussian 
centers in some finite spatial volume what is equivalent to integration over position of one center 
over this volume plus integration over relative distance between centers. The second integration 
includes separation distances between two particles which are comparable with the particle (pion) 
core radius. Indeed, in the pair rest system the probability P2 to registrate two particles with certain 
momenta k„ and k/, can be written in the following form (it is mostly simplified form, for details 



see [b, 16]) 



P2(k«,ki) « J d\ \%/2(x)fJ dx'>D{x,K), (1.2) 
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where D{x,K) is the "relative distance distribution" of the source expressed through the Wigner 
source function S{x,K) (x is relative 4-coordinate), <I>q/2(x) is the wave function distorted by 
particle-particle interaction, which describes a relative motion of two particles, it plays a role of 
a "weight function" or probability to find two pions separated by the distance x with relative mo- 
mentum q = ka — kj„ K = {ka + kh)/2 is a mean momentum of the pair. The integration in (1.2) 
is taken over the 4-volume of the source. If this volume is small the contribution from the region 
where strong interaction dominates appreciable in comparison to contribution from the Coulomb 
interaction of pions. On the other hand, as we shall see in section II, when the emitting volume is 
large one can neglect this region at all because the contribution from strong interaction is relatively 
small. 



The concept of freeze-out accounts for the fact that nevertheless the particles interact always 
electromagnetically one can neglect this interaction in favour of strong interaction in dense many- 
particle system. Then, it is commonly accepted definition that on the background of expansion of 
such a system the beginning of "freeze-out" starts just after a last strong rescattering and in the 
further evolution of the system just electromagnetic interaction left. Natural question can arise in 
this connection - is there contradiction (double counting) with respect to freeze-out definition if one 
includes strong repulsion of pions into account of FSI? To this end let us note, if one integrates over 
the relative distance in (1.2), even from the most naive classical point of view it is obvious, that the 
volume of the pion mean radius {rj^} should be excluded from integration to prevent appearance 
of two pions in the same domain (something like spatial Pauli principal). On quantum mechanical 
level this can be taking into account by the relevant pion-pion repulsive potential. Due to the strong 
interaction the two pions emitted from a small spatial volume, which is comparable with pion 
mean radius (rj^) 0.66 fm, obtain always a repulsive kick which then appreciably reflected in the 

correlation function, as we shall see further. This means that for small volumes (/? < 3 4 fm) 

we cannot sharply separate the electromagnetic and strong interactions and must include the last 
pion-pion strong rescattering in the post freeze-out consideration and by this, rigorously speaking, 
we should reformulate the definition of the "freeze-out". Actually, it is reasonable because the 
registered pions, in some sense, have a memory about their last strong interaction. Indeed, the 
last pion-pion strong rescattering should be accounting separately or extracted from overall strong 
production amplitude because the process of further averaging through multiple rescattering or 
many -particle chaotization is no more going on after it and the last mutual rescattering of registrated 
particles gives a contribution to the post emission relatve dynamics of these particles. It means 
that we separate conventionally the overall creation process on three stages (parts): 1) fireball 
production, 2) strong rescattering of two registered particles on one another, 3) particle-particle 
electromagnetic interaction. This consideration is valid due to inclusive character of two-particle 
spectrum (last strong interactions with other particles are integrated out). Note, that above logics 
deals with quantum mechanical description where one accounts for effects using the probability 
amplitude and two-particle wave function which is spread in space, Othus we can speak about the 
last strong interaction which gives contribution to a certain probability even in the case when the 
centers of wave packets are further from one another than {rji). 
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2. Final state interactions 

Two-particle probability to registrate particles with certain momenta and k;, in the smooth- 
ness approximation was calculated as ^ 

P2{Va,Vb)= j d\d^yS{x+l,Pa)S{x-l,Pb) 0(/)|(/'q/2(y-Vfo/)f + (2.1) 
+ 0(-/)|0q/2(y-Vy)|' ± j ^/"xA5(^+|,^)5(x-f,/:)Cq/2(y-v/)0q/2(y-v/), 

where S{X,K) is the single particle Wigner density of the source 

S{X,K) = j d^xe''^' £ pyy Yr {X + f ) V*y {X - f ) (2.2) 
yY 

with as the density matrix of the source which provides averaging over full set of singl-particle 
quantum numbers 7 of the particle wave function xj/y on the freeze-out hypersurface. We introduce 
pair mean momentum K = (p^ + pfo)/2 and relative momentaum of the registrated particles q = 
Pel —Pb- Using the hermiticity of the density matrix one easily shows that S{X,K) real. Here we 
defined also the three velocities 

v=— , Va = — , Vfo = — (2.3) 

t-K tLK 

associated with the observed particle momenta p^, p^, and their average K. The energy Ek is 
defined as Ek = + and in the pair center off mass system it reduces to particle mass 
Ek = tn. 



The functions in Eq. ( |2.2| ) are the eigenstates of the stationary Schrodinger equation 

^(r)<?'q/2(r)=^rel0q/2(r) (2-4) 

(where £"^1 = q^/2/i and r is the relative coordinate) with asymptotic boundary conditions 

lim 0q/2(r)=^'ii^ (2.5) 

where hamiltonian ^(r) reads 

H{r) = -^yl + V{r) (2.6) 

with reduced mass pL = /m/2 for identical particles and potential which describes the interaction of 
the registrated particles after freeze-out. 

In our further evaluations we will not concentrate on the dependence of the two-particle prob- 
ability PiiPcnPb) on particle velocities putting them approximately equal to the velocity associated 
with mean momentum of the pair, i.e. Vq f« v and v/, v. Then, in pair center off mass system, 
where v = 0, Eq. (2.2) reduces to the expression 

^2(q) = j d'^Xad'^XhS{Xa,Pa) S{xb,Pb) | 0q/2 (Xq - X/,) | ^ 

± J d^Xa d'^Xh S{Xa,K)S{xb, K) */2 (Xfo - Tia ) 0q/2 (Xq - Xfo) , (2.7) 
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where the first term on the r.h.s. of this equation expresses the FSI of the nonidentical particles thus 
it does not factorizes into the product of the two single-particle probabilities. It has transparent 
physical meaning: two single-particle probabilities to find particles in the time-space points Xa and 
Xh with certain momenta pa and p/„ which are expressed by S{x,p), is weighted by the probability 
1 0q/2 (Xfl — x/j) I to find these particles with relative distance — Xh and relative momentum q. The 
correlation function is the ratio of the obtained probability /'2(<l) Eq. (2^) to the product of two 
single-particle probabilities, it reads 



C(q) 



^2(q) 



d'^Xa S{Xa,Pa) / d'^^b S {x}, , Pb ) 



(2.8) 



where 4-vectors pa = {q^/4m, q/2) and pb = (q^/4m, — q/2). 

In the noninteracting limit 0q/2(x) — > exp(/q-x/2) one recovers known expression for the 
correlation function 



C(q,K) = l±. 





1 d'^xe''*''S{x,K) 


2 


J d^XaS{Xa,Pa) J 


( d^XhS{xb,Pb) 



(2.9) 



2.1 Application. Coulomb plus strong final state interaction 



To make evaluations of the correlation function (2.8) we take the source function in the form 



S{x,p) 



(2.10) 



Then, in the pair c.m.s (K = 0) the expression for two-particle probability ( |2.7D reduces to the form 



^2(q) 



jd\e-^'l<\<^^„{rf 

d\e-^'-l<r^,^{-r)^^,^{r) 



(2.11) 



where time dependence is integrated out. In the present simple model of S{x,p) ( 2.10| ) the single- 
particle probability to registrate particle with definite momentum k reduces to the pure Boltzmann 
exponent, i.e. = [47r/M^rA'i (m/rf)]^' exp [— G)(k)/rf], where time and spatial dependencies 

are integrated out. Thus, we can write now the correlation function ( |1.1[ ) in the form 

^ '' "*'"'"°^l^--^'-^l' + '/'q*/2(-r)</'q/2(r) 



C(q) 



7r3/2 (2/?o)3 



d\e-'-"l"'-^ 



</'q/2(r)r±0: 



(2.12) 



where the Boltzmann exponents canceled because 2K^ = ft)(k„) + ft)(kfo) and we adopt that the 
Boltzmann factors have the same freeze-out temperature Tf for the two-particle and single-particle 
spectra (the single-particle spectra are taken usually from different collision event than the two- 
particle spectrum). 

The numerical evaluations of the distorted wave function 0q/2(r) is provided by solution 
of the Schrodinger equation with the relevant potential which reflects two-pion interaction. We 
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solved this equation by partial expansion of the wave function for the Coulomb potential en- 
ergy Vcoui('') = cc/r and for the potential energy which accounts for Coulomb interaction plus 
strong interaction of the frozen pion pair Veff{r) = Vcoui('') + Ktr(''), where "strong potential" 
reads Vstrl?") = Vbexp {—mpr)/{mpr) with Vq = 2.6 GeV, nip = 770 MeV. Potential of the strong 
repulsion [ |l6| ] was chosen to match the behavior of the pion-pion scattering phase shifts. In any 
case the use of this strong potential can be considered as a model of the short-range repulsion 
which is possessed by pions. The correlations functions which correspond to these potentials are 
plotted in Figs. 1-4 for different samples of the freeze-out fireball radii: Rq = 2, 4.5, 7.1, 10 fm. 
It is evidently seen in Figures that the contribution of the strong potential is appreciable for small 
radii emission sources, i.e /?o = 2 fm and Rq = 4.5 fm. For the source size which is of the ^°^Pb 
radius Rq = 1 Aim (see Fig. 3) the contribution which can be put in correspondence to the strong 
potential is comparatively small. And for the source of larger radius, for instance Rq = 10 fm, the 
influence of the strong FSI even negligible (see Fig. 4). 

Remind, that in the absence of the two-particle FSI the source function ( 2.10| ) results in the 



correlation function (K = 0) C(q) = 1 ± exp (— q^/^q). For completeness of comparison we plot 
in Figures this correlation function, corrected by the Gamov factor G(|q|), which for boson-boson 
correlations takes form 

C(q) = G{\q\)(\+e-'^-''l), (2.13) 



where 

G(|q|)=|0q/2(r = 0)1^=^1^ (2.14) 

with T] = amji/\(\\. It is seen in Figures 2, 3 and 4 that the finite size of the emission source softens 
the manifestation of the FSI and the 'Gamov factor' tends to overestimate the FSI effects for the 
source of big size (/?o > 4 fm). 

As it seen in Figs. 1, 2 additional repulsive potential (hadron hard core) essentially suppressed 
the correlation function at small relative momentum. Actually, the physical reasons for this are 

transparent: two strongly interacting particles emitted from the volume of /?o ~ 2 3 fm, when 

their "own mean radius", reflected by repulsive strong potential Vstr, is about (r) 0.66 fm, should 
sufficiently "feel" one another through mutual repulsion. On the other hand, when emision zone 
is much larger than the particle mean radius (r), for instance fireball radius is Rq = 10 fm, the 
contribution of the short-range repulsion is negligible as it seen in Fig. 4. As a matter of fact, the 
last can be regarded as undoubt indication that correlations of the particles which are emitted from 
the large separate distances, starting for instance from R> 4 fm, are not practically influenced 
by hard strong repulsion, or in other words, for the particles which are emitted from perepherical 
regions the contribution from short-range interactions is washed out and their correlations are not 
suppressed by strong interaction. 

It is very interesting to put this conclusion in connection with the so called concept of the 
"dominant contribution to the correlation function of the regions of spatial homogeneity". This 
notion states that on the scale of the collective (hydrodynamical) velocities of two local domains 
it is difficult (small probability) to find two particles with approximately equal momenta. (For 
instance, the most noticeable example is two local parts of the fireball which move in the opposite 
directions to one another.) Thus, the contribution to the correlation function from the regions which 
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Figure 1: Correlation function versus relative pion momentum, Rq = 2.0 fm. Dotted (far-between) curve 
is the correlation function without FSI. Fat solid curve corresponds to the two-particle Coulomb potential 
Ucouli'') = Oi/r. Dotted curve corresponds to the sum of the Coulomb and strong two-particle potentials 
f/ (r) = a/r -h Vb exp {—mpr)/ {ntp r) . The bottom solid curve is the correlation function corrected just by the 
Gamov factor. 

spatially stand far from one another is suppressed due to collective motion. Actually, in the present 
simple model we do not consider expansion of the source. 

On the other hand, as we have just proved a short-range repulsion of the strongly interacting 
particles emitted from a small domain (/? < 4 fm) decreases as well the probabihty to find two 
particles with a small relative momentum, for example if the volume of homogeneity is of the size 
/? PS 2 — 4 fm, then the correlation function of pions from this region, as we see from Figs. 1, 2, will 
be suppressed by strong repulsion of these particles. Hence, both effects work in the same direction, 
decreasing the correlations of the particles emitted from the small regions (hard core repulsion) and 
decreasing the correlations of the particles emitted from the far separated regions (collective mo- 
tion). If these effects have the same scale then they approximately uniformly decrease correlations 
of two particles which are emitted from all regions of the fireball. The last statement results in con- 
clusion that two-particle correlations does not occur preferably from the small regions where the 
spatial separation of the particles is small (hence, they have the same collective velocity), but the 
particles which are emitted from two regions which are separated by large distance contribute also 
appreciably to the correlation function. Thus, the correlation function shows us the real size of the 
fireball at freeze-out, but not just a size of the small region of homogeneity. If the region of spatial 
homogeneity has the radius /? > 7 fm then the strong repulsion can be neglected. So, the strong FSI 
effectively extends the size of the region from which the particles are allowed for correlations. If 
we unite two cases where the radius of homogeneity is small /?homo < 4 fm and large /?homo > 7 fm 
(hence, it is comparable with the size of the fireball) then we come to the conclusion that the radius 
which is extracted from the correlation function is approximately the radius of the fireball. 

It is interesting to point out that for the sources of the size /?o ~ 3 fm correction of the cor- 
relation function just by the Gamov factor effectively gives a right correction, Uterally it becomes 
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Figure 2: The same as in Fig. 1, but for Rq = 4.5 fm. 
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Figure 3: The same as in Fig. 1, but for Rq = 1 A fm. 

equal to the account of a finite size of the source exploiting potential Veffir) which includes strong 
FSI. It means, that if one starts from the Gamov factor correction then account for finite size of the 
source softens Gamov's correction. If we include strong interaction (in addition to Coulomb one) 
then for the source of the small size account for finite size of the source makes correction more 
hard (in comparison to Gamov factor correction). But in the intermediate point, which is /?o ~ 3 fm 
these two deviations from the Gamov corrections cancel one another 



3. Summary and conclusions 

We examined pion-pion final state interaction when the two-particle potential energy is not just 
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Figure 4: The same as in Fig. 1, but for Rq = 10 fm. 



the Coulomb one. Distortion of the Coulomb potential is taken in the additive form, what is due to 
account for particle-particle strong interaction. To make our investigation transparent as much as 
possible we take: 

a) for the description of the fireball (particle emitting source), the Wigner source function S{x,K) in 
the Gaussian form; b) small relative momentum approximation what results in equlity of two parti- 
cle velocities and pair c.m.s. velocity, then expression for the two-particle probability is sufficiently 
simplified, see Eq. (2.7). 

We argue that for inclusive two-particle spectrum one should explicitly account for strong in- 
teraction of two registered particles because their mutual last rescattering essentially determines 
dynamics of their post-freeze-out relative motion. Strong two-particle potential was chosen as 
screened Yukawa potential and reflects short-core repulsion of pions. Then the wave function de- 
scribing the relative motion of pions was founded as numerical solution of the Schrodinger equation 
where potential is taken as sum of the Coulomb and Yukawa potentials. We obtained that corre- 
lation function is sensitive to the presence of short range two-particle interaction when emission 
volume is of the radius Ro <4 fm, see Figs. 1 and 2. Indeed, two strongly interacting pions emitted 
from the separate distance (separation of the centers of the wave packets), which is comparable 
to the pion "mean radius" (rji) 0.66 fm, obtain a mutual repulsive kick which is then reflected 
as decreasing of the two-particle probability to find particles with "initial production" relative mo- 
mentum, i.e. the "strong kick" shifts along x-axis this particular probability value to higher values 
of the relative momentum. Consequently, this will result in suppression of the correlation function 
what is most noticeable for small relative momentum. (Obviously, the effect which is due to the 
strong repulsion is additional to that one which is caused by the Coulomb repulsion.) 

On the other hand the presence of strong particle-particle repulsion almost unnoticeable in 
the correlation function when one considers the sources of big size, Rq >^ fm, see Figs. 3 and 
4. This just confirms that the hard core repulsion of the registered particles appreciably influences 
the correlation function when the particles are emitted from the region of a small size. If the 
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region of homogeneity is also small /?homo < 4 fm then the correlations originated by pions emitted 
from this region are suppressed. This suppression can be of the same strength as suppression of 
the correlations of pions emitted from the domains which are beyond the radius of homogeneity 
(the latter suppression is due to collective motion of expanding system) what results in uniform 
suppression of correlations of two pions emitted from Then, the radius Rq which is measured by 
the correlation function can be larger than the size of homogeneity region and tends to cover the 
overall size of the fireball. In the opposite case, when the size of the region of homogeneity is 
larger /?homo > 7 fm, the correlation function is not influenced practically by the strong interaction 
of two registered pions. Hence, the radius Rq which is extracted from the correlation function in 
this case is about /?homo and Rq really reflects the size of homogeneity volume. 

Summarizing all above we can say: 
The strong final state interactions can be neglected when the volume of spatial homogeneity in the 
fireball are of the size /?homo > 7 fm. Then, the radius Rq which is extracted from the pion-pion 
correlation function measures the size of the regions of homogeneity in the emitting source and 
coincide with /?homo- When the radius of homogeneity volume is small /?homo < 4 fm the strong 
final state interactions are noticeable and their presence is reflected in the correlation function in 
the way that radius Rq (extracted from the correlation function) is larger than /?homo and tends to be 
the overall fireball radius. 
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